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ABSTRACT: We introduce a model to describe tumor development. The model takes into account the 
development of their own cells, the effects of chemotherapeutic drugs on them, as well as their 
interaction with surrounding cells. We also introduce an analytical approach for the analysis of the 
introduced model. The approach gives the possibility to take into account changes in the conditions of 
the processes considered: changes in the properties of tissues in space and time, as well as the 
nonlinearity of the development of tumors. We analyzed the possibility of accelerating and 
decelerating the development of the tumor depending on the parameters of the considered process. 
The possibility of changing their rates is being considered. 
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INTRODUCTION 


The development of cancerous tumors is associated with various chemical, genetic, physiological, 
and mechanical factors that occur both at the subcellular and cellular levels, as well as at the levels of 
tissues and organs. In recent decades, one can find progress in identifying and explaining the processes 
that occur during the development of cancer, as well as developing methods and tools for early 
diagnosis and treatment of the disease [1-7]. The development of biotechnologies and medicine has 
allowed us to lay significant groundwork for the treatment of the disease in question. A significant part 
of the study of tumor development is its prognosis based on experimental data and mathematical 
models. We introduce a model of tumor development. The model gives them the possibility of taking 
into account the development of their own cells, the effects of chemotherapeutic drugs on them, as 
well as their interaction with surrounding cells. We also introduce an analytical approach for the 
analysis of the introduced model. The model gives the possibility of taking into account changes in the 
conditions of the considered processes. We also consider the possibility of changing their flow rate. 
Method of solution 


In this section, we consider a model to estimate the development of tumors, their own growth, the 
influence of chemotherapeutic drugs on the above cells, as well as their interaction with surrounding 
cells. We determine the required spatio-temporal distribution of the number of the considered cells as 
the solution of the second Fick's law in the following form: 
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where N (x,y,z,t) is the quantity of tumor cells; terms with functions Rj(x,y,z,t) describes the rate of 
division of tumor cells; terms with function M (x,,z,t) describes the rate of natural mortality of tumor 
cells; the term with the functions Aj(x,y,z,t) describes the rate of death of tumor cells due to their 
bioavailability to an active chemotherapeutic drug; S (x,y,z,t) is the rate of death of tumor cells due to 
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competition with surrounding cells; D (x,y,z,t) is the cell diffusion coefficient; x, y and z are the spatial 
coordinates; t is the time; T is the temperature. The initial distribution of tumor cells and the boundary 
conditions for equation (1) are represented in the following form 
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N (xy,z,0)=f (x,y,z). 

We determine the solutions of equation (1) using the method of averaging functional corrections [8- 
11]. In the framework of the method, we will substitute in the right parts of equations (1) the not yet 
known average value of the desired number of cells a instead of the function N (x,,z,t), i.e. we use the 
following replacement N (x,y, z,t) >a1. Then we obtain the following equation to determine the first- 
order approximation of the number of tumor cells in the following form: 
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Integration of the right and the left sides of equation (3) on time gives a possibility to obtain a ratio 
for the first-order approximation of the number of the considered cells Ni(x,y,z,t) in the final form 
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The not yet known average value a was calculated by using the following standard relation [8-10] 
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Substitution of relation (4) into relation (5) gives a possibility to obtain the following equation to 
determine the above average value 
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Solution of equation (6) could be written as 
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The second-order approximation of the az + N; (x,y,z,t) in the right side of equation (1) 
required quantity of tumor cells in the [8-11]. The substitution gives a possibility to 
framework of the method of averaging function obtain the following equation to calculate the 

corrections could be obtained by using the considered approximation 
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Integration of right and left sides of equation (7) on time gives a possibility to obtain the following 


relation to calculate the second-order approximation of the considered cells N2(x,y,z,t) in the following 
final form: 
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Average value of the second-order approximation of the considered quantity of sells a2 could be 
calculated by using the following standard relation [8-11] 
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Substitution of relations (4) and (8) into relation (9) gives a possibility to obtain the following 
equation to determine the not yet known average value az 
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Solution of equation (10) could be presented in the following form 
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The spatial-temporal distribution of the 
quantity of tumor cells was analyzed 
analytically by using the second-order 
approximation in the framework of the method 
of averaging function corrections. The 
approximation is usually a good enough 
approximation to make a qualitative analysis 
and obtain some quantitative results. All the 
obtained results have been checked by 
comparison with the results of numerical 


simulations. 


DISCUSSION 

In this section, we analyzed the changes in 
the distribution of the quantity of tumor cells 
in an organism over space and time. Figure 1 
shows typical dependences of the considered 


1.00 


0.75 


N(x.) 


0.25 


-2.50 


Af f fa (x,y, z,t) dzdydxdt 
000 


E-ISSN: 2349 5359; P-ISSN: 2454-9967 


L, Ly 


“x 


a (x, y,z,t)R, (x, y, zt) C(x, y, zt) x 
000 
GE 


L 


“yy 


o-af ffm (x, y,z,t)d zd ydxdt+L,L,x 
000 


& (2) 
=f NY {4.6 x, y,Z,t)R, (x, y,z,t) N; (x, y,z,t)dzd ydx eae 


L, 
[f4, (x, y, z,t)R, (x, y,z,t) [N, (x, y,z,t)+1]dzd ydxdt [L,L,L.O+ 
00 


(11) 


quantity on coordinates at constant values of 
parameters. Depending on the values of the 
parameters in equation (1), the tumor may 
grow or shrink. The rate of change in the size of 
the tumor could also be changed. The 
effectiveness of the chemotherapy drug 
increases when combined with an external 
low-intensity physical factor that increases the 
biological activity of both tumor cells and 
tumor tissue. External physical factors become 
with their spatial inhomogeneity. Decreasing of 
size of the neoplasm during treatment in the 
framework of the considered was observed 
when the intensity of the external field is above 
a certain threshold level. The similar 
conclusion was qualitatively obtained in the 
following references [12,13]. 
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Fig. 1. Typical normalized dependences of quantity of tumor cells on co-ordinate 


CONCLUSION 

We introduced a model for the description of 
the development of tumors. The model gives 
the possibility to take into account the 
development of the tumor's own cells, the 
effects of chemotherapeutic drugs on them, as 
well as their interaction with surrounding cells. 


We also introduced an analytical approach for 
the analysis of the considered model. The 
approach gives the possibility of taking into 
account changes in the conditions of the 
considered processes. We also consider the 
possibility of changing their rates. 
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